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Q ■ Abstract 

In this paper we introduce the exponentiated Weibull power series (EWPS) class of distributions 

■ which is obtained by compounding exponentiated Weibull and power series distributions, where 
I the compounding procedure follows same way that was previously carried out by Roman et al. 

(2010) and Cancho et al. (2011) in introducing the complementary exponential-geometric (CEG) 
i and the two-parameter Poisson-exponential (PE) lifetime distributions, respectively. This dis- 

^ tribution contains several lifetime models such as: exponentiated weibull-geometric (EWG), 

exponentiated weibull-binomial (EWB), exponentiated weibull-poisson (EWP), exponentiated 
' weibull-logarithmic (EWL) distributions as a special case. 

, The hazard rate function of the EWPS distribution can be increasing, decreasing, bathtub- 

I shaped and unimodal failure rate among others. We obtain several properties of the EWPS 

. distribution such as its probability density function, its reliability and failure rate functions, 

quantiles and moments. The maximum likelihood estimation procedure via a EM-algorithm is 

■ presented in this paper. Sub-models of the EWPS distribution are studied in details. In the 
I end. Applications to two real data sets are given to show the flexibility and potentiality of the 

EWPS distribution. 

Keywords: EM algorithm. Exponentiated Weibull distribution. Maximum likelihood 
estimation. Power series distributions. 
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1. Introduction 

The Weibull and exponentiated Weibull (EW) distributions in spite of their simplicity in 
solving many problems in lifetime and reliability studies, do not provide a reasonable parametric 
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fit to some practical applications. 

Recently, attempts have been made to define new families of probability distributions that 
extend well-known families of distributions and at the same time provide great flexibility in 
modeling data in practice. One such class of distributions generated by compounding the 
well-known lifetime distributions such as exponential, Weibull, generalized exponential, expo- 
nentiated Weibull and etc with some discrete distributions such as binomial, geometric, zero- 
truncated Poisson, logarithmic and the power series distributions in general. The non-negative 
random variable Y denoting the lifetime of such a system is defined by y = mini<j<jvXj or 
Y = maxi<j<ArXj , where the distribution of Xi belongs to one of the lifetime distributions and 
the random variable N can have some discrete distributions, mentioned above. 

This new class of distributions has been received considerable attention over the last years. 
The exponential geometric (EG), exponential Poisson (EP), exponential logarithmic (EL), ex- 
ponential power series (EPS), Weibull geometric (WG), Weibull power series (WPS), exponen- 
tiated exponential- Poisson (EEP), complementary exponential geometric (CEG), two-parameter 
Poisson-exponential, generalized exponential power series (GEPS), exponentiated Weibull-Poisson 
(EWP) and generalized inverse Weibull-Poisson (GIWP) distributions were introduced and stud- 

Barreto-Souza et al. [Tj, Morals and Barreto-Souza et al. [23(, Barreto-Souza and Cribari- 
Neto Is]; Louzada-Neto et al. 1^, Cancho et al. [lo| . Mahmoudi and Jafari [3], Mahmoudi 
and Sepahdar [20| and Mahmoudi and Torki [21]. 

In this paper we introduce the exponentiated Weibull power series (EWPS) class of distribu- 
tions which is obtained by compounding exponentiated Weibull and power series distributions, 
where the compounding procedure follows same way that was previously carried out by Roman 
et al. (2010) and Cancho et al. (2011) in introducing the complementary exponential-geometric 
(CEG) and the two-parameter Poisson-exponential (PE) lifetime distributions, respectively. 
This distribution contains several lifetime models such as: exponentiated weibull-geometric 
(EWG), exponentiated weibull-binomial (EWB), exponentiated weibull-poisson (EWP), expo- 
nentiated weibull-logarithmic (EWL) distributions as a special case. 
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2. Exponentiated Weibull distribution: A brief review 

Mudholkar and Srivastava {24] introduced the EW family as extension of the Weibull family, 
which contains distributions with bathtub-shaped and unimodal failure rates besides a broader 
class of monotone failure rates. One can see Mudholkar et al. [2511, Mudholkar and Huston 
Gupta and Kundu [l^, Nassar and Eissa 



26] 



28( 1 and Choudhury [13] for applications of the EW 
distribution in reliability and survival studies. 

The random variable X has an EW distribution if its cumulative distribution function (cdf) 
takes the form 

Gx{x) = (l-e-^^^y'y, x>0, (1) 

where 7 > 0, a > and /3 > 0, which is denoted by EW{a, /3, 7). The corresponding probability 
density function (pdf) is 

a-l 



gx 



(x) = a7/3^x^-ie-(^^')" (l 



(2) 



The survival and hazard rate functions of the EW distribution are 

5(x) = l-(l-e-(^-)y, 



and 



h{x) = a-f^^x^-^e-^'^''^'' (l - e'^^^^ 



0-1 



i_fi_e-(/^-r 



respectively. The kth. moment about zero of the EW distribution is given by 



EiX'') = a(3-''T ( - + 1 ) 



7 



j=0 



vr-ij(, + i)-(^i). 



(3) 



Note that for positive integer values of a, the index j in previous sum stops at a — 1, and the 
above expression takes the closed form 

'k 



E{X') = ar'r ( - + 1 ) Akij) 



where 



a-l 



-(^+1) 



A; = 1,2,3,' 



(4) 



(5) 



in which r(^ + 1) denotes the gamma function (see, Nassar and Eissa (2003) for more detail). 
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Table 1: Useful quantities of some power series distributions. 



Distribution 




C{9) 


c'{e) 


C"{0) 




S 


Poisson 


n!-i 








log ((? + !) 


oo 


Logarithmic 


71-1 


-log (1-0) 








1 


Geometric 


1 




(1-9)-' 


2(l-6l)-3 


9{l + 9)-^ 


1 


Binomial 


o 


(61 + 1)"-! 




m{7n—l) 


{9-1)^-1 


oo 



3. The class of EWPS distribution 

Consider the random variable X having the EW distribution where its cdf and pdf are given 
in ([I]) and 

Given A^, let Xi,--- ,Xjv be independent and identically distributed (iid) random variables 
from EW distribution. Let the random variable is distributed according to the power series 
distribution with pdf 

where an > depends only on n, C{0) = Ti^^iUnO'^, £ (0, s) is such that C{6) is finite. For 
more details on the power series class of distributions, see Noack (1950). Table 1 shows useful 
quantities of some power series distributions (truncated at zero) such as poisson, logarithmic, 
geometric and binomial (with m being the number of replicas) distributions. 
Let Y = max(Xi, • • • ,X]\[), then the conditional cdf of Y\N = n is given by 



GY\N=n{y) = {G{y)r = (l - e-(^^)^)"" , (6) 



which is the EW distribution with parameters na, (3, 7, and denoted by EW(na, /3, 7). The ex- 
ponentiated Weibull power series (EWPS) distribution, denoted by EWPS (a, (3, 7, 9), is defined 
by the marginal cdf of Y, i.e.. 



Remark 1. Let Y = min(Xi, • • • , Xjv), then the cdf of Y is given by 



7 \ ce 



FY{y) = i ^ y>o. (8) 

If a = 1, then the cdf of Y is Fyyu) = 1 ^ ^ ^^^^ which is called Weibull Power Series 

distributions (Morals and Barreto-Souza, 2011) and this family includes the life time distribution 



presented by Barreto-Souza et al. (2010a), Barreto-Souza et al. (2010b). which Xj's has the 
exponentiated Weibull distribution is obtained. The EG distribution (Adamidis and Loukas, 
1998) is obtained by taking C{e) = 9(1 - 9)-^ with 6^ G (0, 1) and a = 1, 7 = 1 in ([8]). Moreover, 
for a = 1,7 = 1, we obtain the EP distribution (Kus, 2007) and the EL distribution (Tahmasbi 
and Rezaei, 2008) by taking C{9) = e^-l,90, and C{9) = - log (1 - 6*), 6* G (0, 1), respectively. 
The WG distribution (Barreto-Souza et al. (2010a), Barreto-Souza et al. (2010b) ) is obtained 
by taking C{9) = 9{1 - 9)'^ with 9 £ (0, 1) and a = 1 in dS]). The EWG distribution is obtained 
by considering C{9) = 9 {I - 9)'^ with 9 G (0, 1) in ([8]). 



The pdf of the EWPS distribution is given by 

Ci9) 



.^.,C'(9(l-e-(Py)''r) 
fviy) = 0a7/3V-'e-(/^J^)^ (l - 6-^^-^) ^, (9) 



where a, /3, 7 > and 9 G (0, s). 

The survival function and hazard rate function of the EWPS distribution are given, respectively, 

by 

C(9{l-e-(^y^''r) 
S{y) = l-^ C{9) ' 

and 

C'(9(l-e-(^yy'r) 

h(v) = 9a'yB'^v'^"^e~'^'^y'i^' ( 1 - e'^'^^^" I — fill 

n^yj va^p y yi e ^ C{9) - C {9 {l - e-i^y)-")")' ^ ' 

Consider C{9) = 9 + 9"^^. If /? = 1 and 9 = 1, the plots of this density and its hazard rate 
function, for a = 0.5, 7 = 1,0 = 1,7 = 0.5, q = 2, 7 = 3, and a = 2, 7 = 1 are given in Fig 1. 

Proposition 1. The limiting distribution of EWPS{a, (3,^,9) when — )■ 0^ is 

hm ^F(„\ -lim ^ C{eG(y)) _ ,. E^T^i "ne"(G(y))" 

- 11,^ «c(G(y))-+E;r=.+iane"--(G(s/))" 



= {G{y)Y = {l-e-^Py)''y'' 

which is a EW distribution with parameters ca, 7 and /3, where c = min{n G N : Un 0}. 

Proposition 2. The densities of EWPS class can be expressed as infinite linear combination 
of density of order distribution. We know that 



C'{9) = T^^^^nan^ 



jn— 1 




Figure 1: Plots of pdf and hazard rate function of EWPS where C(6) = 6 + 6' 



Therefore, 

71=1 

Using the EW density given before, we obtain 

fEWPsiy; a, /3, 7, 0) = 6^ E -^r^fEwiy; na, 13, 7). (12) 

n=l 

Various mathematical properties (cdf, moments, percentiles, moment generating function, fac- 
torial moments, among others) of the EWPS distribution for \9\ < 1 can be obtained from Eq. 
([12]) and the corresponding properties of the EW distribution. 

Proposition 3. The density of EWPS distribution can be expressed as infinite linear com- 
bination of density of the biggest order statistic of Xi,--- where Xi ~ £'VF(a,/3,7) for 
i = 1,2, ■ ■ ■ ,n. we have 

00 00 
fEWPsiy) = Y.{G{y)rP{N = n) = J29x,Jy)P{N = n), 

n=l n=l 

in which gx^^^ {y) is the pdf of = max(Xi, • • • , X„). 
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4. Quantiles and moments of the EWPS distribution 

The qth quantile of the EWPS distribution is given by 

Where G~^{y) = ^ ^— In ^1 — ^ , and C^^(.) is the inverse function of C(.). The qth 
quantile of the EWPS distribution is used for data generation from the EWPS distribution. In 
particular, the median of the EWPS distribution is given by 

' C-^{0.5C{e)) \ 



-1 



yo.5 = G ^ Q J 

Suppose that Y ~ EWPS{a,l3,j,9), and = max(Xi, • • • where Xi ~ EW{a, (3,^) 

for i = 1,2, ■ ■ ■ ,n, then the kth moment of Y is given by 

E{Y'^) = E{E{Y^\N)) = EZi m = n)E{Y'^\n) = E~ i = n)E{Xf^^) 

= ap-'T U + l) En=i Er=o(-lP?CrT>(j' + ^r^^^'^- (13) 



For positive integer values of a, the index j in above expression stops at na — 1 

Now we obtain the moment generating function of the EWPS distribution using the Eq. 
([13]), as follow 

My(t) = E.=ol^(n 

= E£o I [(1 - ^)«/3-T (i + 1) En=i Er=o(-iy rTV^""'(^' + ^y^^^'\ 

_ v^oo v^oo v^oo i-iy /, /oy a„e" (na-l\ "-"^(^+1) 

— 2^1=0 2^n=i 2^j=o ii v-IP) c{e)\ j ) ^ \^^) 

According to the Eq. (I13p . the mean and variance of the EWPS distribution are given, respec- 
tively, by 



i^(y) = ar'r(- + l)>^>J-l)^^{ )n(J + lr--^^ (15) 

and 

Var{Y) = ar'T + l) E E(-l)^'^ ^ ') n{j + 1)"^^^) - E\Y), (16) 

' 71=1 J=0 \ I \ J / 

Where -E(y) is given in Eq. (llSp 
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5. Renyi and Shannon entropies 

Entropy has been used in various situations in science and engineering. The entropy of a 
random y is a measure of variation of the uncertainty. For a random variable with the pdf /, 
the Renyi entropy is defined by lR{r) = \og{ f^f^{y)dy}, for r > and r 7^ 1. For the 
EWPS distribution, the power series expansion gives 



c{e) 



^.,,^_, [C'(e(l-e-^M-'r)Y 

[c(e)Y 



Applying the Equation [Yll^o'^iU^) ^ = Y,i=o^i,j^ > v,x.v..v. ...^ ^^^...^.^..^^ 



°° Ci.iU*, where the coefficients c,;., for i = 1,2, 
can be easily obtained from the recurrence relation Cij = {iwo)~^ Ylm=iU^ — i + m)wmCi-m,j , 
white CO J = wo^ for [C'{9{1 - e'^^?')^)")] '^and series expansion for (l - (1 - e-(/'^)^)")"^''"^^^ 
gives 



c(e) 



But setting u = (j + r){(3y)^, gives 

Io^ny)dy 

Substituting from ()17p . we obtain 
1 



ad 



C(9) 



U+r) 



00 (~iye^fa{r+i)- 



lR{r) 



1 — r 



■log < 







< 


[c{e)\ 



7 



^ 00 CXD 



(-1)^(9* /a(r + i) - r 



i=o j=o {j + r)' 

The Shannon entropy which is defined by E[— log(/(y))], is derived from limr_s.i iRir). 



(17) 



(18) 



6. Moments of order statistics 

Order statistics make their appearance in many areas of statistical theory and practice. Let 
the random variable Yr;n be the rth order statistic (Yi-n ^ Y2:n < • • • < Yn-.n) in a sample of size 
n from the EWPS distribution. The pdf of Yr-.n for r = 1, • • • , n, is given by 

fr.niy) = ^(^^ ^ ! , ^ 1) mFiyy-' [1 - F{y)r-^, y>0. (19) 
where F{y) and f(y) are given in ([7]) and ([9]). Substituting from ([7]) and ([9]) into (fT9|) gives 

{c{e{i - e-^^y^''y)Y-\c{9) - c{e{i - e-(/5?^)^)"))"-^ (20) 
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Also the cdf of is given by 

Fr:n{y) = EL. (fc) [^(y)]' " ^(y)]""' 

— l^k=r [kJ c(ey' ■ ^'^^> 

Expression for the rth moment of the order statistics (Yi-n ^ ^2:n < • • • < Yn-n), with a 
cdf in the form (j2ip . are obtained by using a result due to Barakat and Abdelkader (2004) and 
becomes 



= ^Er=n-.+i Er=o ^^i^^S^O r y'-'Cm - e-i^yVTYdy. (22) 



7. Residual life function of the EWPS distribution 

Given that a component survives up to time t > 0, the residual life is the period beyond t 
until the time of failure and defined by the conditional random variable X\X > t. In reliability, 
it is well known that the mean residual life function and ratio of two consecutive moments of 
residual life determine the distribution uniquely (Gupta and Gupta, 1983). Therefore, we obtain 
the rth order moment of the residual life via the general formula 

mr{t) = E[{Y- tr\Y >t\ = ^^j^{y- tYf{y)dy, 

where S{t) = 1 — F{t), is the survival function. 

Applying series expansion (9), the binomial expansion to {y — ty and substituting S{t) given by 
(jlOp into the above formula gives the rth order moment of the residual life of the EWPS as 

^"■^^f ~ C(e)-C(e(l-e~(W^)") l^i=0l^n=ll^j=Q ^,^_^^r_+jp^\ j ) 

x$(^;(i + l)(/3t)^), r>l, (23) 

where <I>(s; t) is the upper incomplete gamma function given by <I>(s; t) = x^~^e~^dx. 
Another important representation for the EWPS is the mean Residual life (MRL) function 
obtain by setting r = 1 in Eq. (I23p . MRL function as well as failure rate (FR) function is 
very important since each of them can be used to determine a unique corresponding life time 
distribution. Life times can exhibit IMRL (increasing MRL) or DMRL (decreasing MRL). MRL 
functions that first decreases (increases) and then increases (decreases) are usually called bathtub 
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(upside-down bathtub) shaped, BMRL (UMRL). The relationship between the behaviors of the 
two functions of a distribution was studied by many authors such as Ghitany (1998), Mi (1995), 
Park (1985), Shanbhag (1970), and Tang et al. (1999). For the EWPS distribution the MRL 
function is given in the following theorem. 

Theorem 1. The MRL function of the EWPS distribution with cdf ^ is 

m{t) = {fii + I{t)-t)/Sit), t>0 (24) 

where I{t) = F{y)dy, S{t) is the survival function given in [W\). and is the mean of the 
EWPS in Eq. (d). 



Proof. For more detail about the proof of this theorem see Nassar and Eissa (2003). 
According to Theorem 1, for the EWPS distribution with f{y) given by Q, we have 

where ^{s]t) is the lower incomplete gamma function given by^'(s;t) = jlx'-^e-''dx. Substi- 
tuting Eqs. ([El), ^ and ^ into ([MD gives the MRL function of the EWPS distribution. 



{-If an 



'mi{t) — C{e)l3^n=l^k=0 s{t) 

8. Reversed residual life function of the EWPS distribution 

Given that a component survives up to time t > 0, the residual life is the period beyond t 
until the time of failure and defined by the conditional random variable X\X > Therefore, we 
obtain the rth order moment of the residual life via the general formula 



mr{t) 



E [{Y - tr\Y >t] = ^^ fiy - tyfiy)dy, 



where F(t), is The exponentiated Weibull power series (EWPS) distribution. 

Applying series expansion (9), the binomial expansion to {t — yY and substituting F{t) given by 

([7|) into the above formula gives the rth order moment of the reversed residual life of the EWPS 

as 

— a v^r v^oo v^oo 'na„6l" (na-l\ 

Mr(t) - 2^i=Q l_,n=i 2^j=o~ ,+1 , , . \ j ) 



xM/(^l + ^;(j + l)(/3t)^j, r>l, (26) 
where ^{s]t) is the upper incomplete gamma function given by $(s;t) = x^~^e~^dx. 
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9. Probability weighted moments 

Probability weighted moments (PWMs) are expectations of certain functions of a random 
variable defined when the ordinary moments of the random variable exist. The probability 
weighted moments method can generally be used for estimating parameters of a distribution 
whose inverse form cannot be expressed explicitly. We calculate the PWMs of the EWPS 
distribution since they can be used to obtain the ordinary moments of the EWPS distribution. 

The PWMs of a random variable X are formally defined by 



oo 



Ts,r = E[X'FiXY] = / x'F{xYfix)dx, (27) 
Jo 

where r and s are positive integers and F{x) and f{x) are the cdf and pdf of the random variable 

X. The following theorem gives the PWMs of the EWPS distribution. 

Theorem 2. The PWMs of the EWPS distribution with cdf ^ and pdf are given by 



Proof. Substituting from ^ and dH]) into ([27]) gives 



oo 




Applying the Equation ( X]£o '"^«'"*) ~ Z]£o^*.i^*' where the coefficients for i = 1,2,... 
can be easily obtained from the recurrence relation Cij = {iwo)~^ X^m=iO'"^ — i + m)wmCi-m,j , 

white Co J = wq^ for [C"(0(1 — e"^''^)'')")] ''and series expansion for (l — (1 — e"^''^^^)'^) 
gives 



a^e^l oo oo oo 

T. ^ 



n=l 1=0 7=0 ^ ^ n 



ic{e)y+ 

Setting u = {j + l){l3xY the above integral reduces to 







oo 

^s+^-i^-ik+iKM-dx = — -^r(i + -) 



and the proof is completed. 

Remark 2. The sth moment of EWPS distribution can be obtained by putting r = in Eq. 
(1281) . Therefore 

Eixn = ail - e)r'ri^ + 1) E E(-i)'=a + 1) ("^'' ~')ik + 1)-^^^'^9\ 
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which is equal with Eq. (jl3p if s is replaced by k. Also, the mean and variance of the EWPS 
distribution can be obtained. 



10. Mean deviations 

The amount of scatter in a population can be measured by the totality of deviations from 
the mean and median. For a random variable X with pdf /(.), cdf F{.), mean /i = E(X) 
and M = Median{X) the mean deviation about the mean and the mean deviation about the 
median, respectively, are defined by 

5i{X)= / \x- ii\f{x)dx = 2iiF{^i)-2ii + 2L{p), 
Jo 



and 



/•oo 

62iX)= / \x-M\f{x)dx = 2MF{M)-M-fi + 2L{M), 
Jo 



where L(^) = f°° xf{x)dx and L{M) = xf{x)dx. 



For the EWPS distribution we have 



n=l j=0 \ I \ J J V -r 7 / 

Theorem 3. The Mean deviations of the EWPS distribution are 

/c(9(l-e-(f'')'')°) -C(9)\ 

= [- cW) j + 

and 

52 = 2L{M) - fi, 

respectively, where fi is the mean of EWPS in Eq. il5\) . L{fi) and L{M) are obtained by substi- 
tuting fi and M = j^ c-i(o^5C(e)) ^ instead of b in Eq. <fIP|). 

11. Bonferoni and Lorenz curves 

The Bonferroni and Lorenz curves and Gini index have many applications not only in eco- 
nomics to study income and poverty, but also in other fields like reliability, medicine and insur- 
ance. The Bonferroni curve Bf[F{x)] is given by 

Bf[F{x)] = —3— ruf{u)du. 
f^F{x) Jo 
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The Bonferroni curve of the EWPS distribution is given by 

x(i + l)-^^+')^((i + l);(j + l)(x/3r). 
Also, the Lorenz curve of EWPS distribution can be obtained via the expression 

x(i + l)-(^i)*((i + l);(i + l)(x^r). 

The scaled total time on test transform of a distribution function F (Pundir et al., 2005) is 
defined by 

1 /•* - 
SF[F{t)] = - F{u)du. 

Jo 

If F{t) denotes the cdf of EWPS distribution then 

smt)\ = - ^Efi-D'^M* ( -■•my) . 

The cumulative total time can be obtained by using formula Cp = Jq SF[F{t)]f{t)dt and the 
Gini index can be derived from the relationship G = 1 — Cp- 

12. Estimation and inference 

In what follows, we discuss the estimation of the parameters for the EWPS distribution. Let 
^1) ^2, • ■ ■ , ^ be a random sample with observed values yi,y2, - ■ ■ ,yn from EWPS distribution 
with parameters a,P,j and 9. Let = {a,P,j,9)'^ be the parameter vector. The total log- 
likelihood function is given by 

In = Iniv, e) = n[loga + log7 + 7log/9 + log^] + (7 - 1) ELi ^^^Vi ' Y^tMiV 

+ («-!) Ei=i log(l - e-^f^y^^') -nlogC{e) + EILi logic (9 (1 - e-(^f^)'')")]. 

The associated score function is given by Un{@) = {dln/da,dln/dp,dln/dj,dln/dX)'^ , where 
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din 
dp 



din 

97 



9k 
961 



-I- ca7P Z^i=i c'(e(i-e-('»!'.) ')") ' 

+ I" - J-J Z^i=l i_e--{fiy^)^ ^ 2-1=1 C'(e(l-e-('5!'i)'')«) ' 



n 1 Y^n 



n (l-e-tfi'»)^)°C"(6'(l-e-('^i^»)^)°') 



— n 



C"(6»(l-e-(^2'<)'^)«) 

The maximum likelihood estimation (MLE) of 9, say 0, is obtained by solving the nonlinear 
system [/„ (G) = 0. The solution of this nonlinear system of equation has not a closed form. For 
interval estimation and hypothesis tests on the model parameters, we require the information 
matrix. The 4x4 observed information matrix is 



4(e) = - 



laa Ia0 ^a'y laO 

lap 1/3(5 Il3"i 1/36 
-^a7 -^/97 -^77 ^-fO 

lao Ipe I^e he 

whose elements are given in Appendix. 

Applying the usual large sample approximation, MLE of 6 i.e. can be treated as being 
approximately N/^{Q, Jn{Q)~^), where J„ (0) = E[In (©)]■ Under conditions that are fulfilled 
for parameters in the interior of the parameter space but not on the boundary, the asymptotic 
distribution of ^(9-6) is Ni{Q, J{Q)~^), where J (9) = lim 

n— ^00 ^ -^ri(9) is the unit in- 
formation matrix. This asymptotic behavior remains valid if J(9) is replaced by the average 
sample information matrix evaluated at 9, say n~^/„(9). The estimated asymptotic multivari- 
ate normal A/4(0,7„(9) ^) distribution of can be used to construct approximate confidence 
intervals for the parameters and for the hazard rate and survival functions. An 100(1 — 7) 
asymptotic confidence interval for each parameter 9r is given by 



ACIr = (0r- - Z2. V Qr + ZxV Z^^, 



where J^'" is the {r, r) diagonal clement of In ^9^ 
1 — 7/2 of the standard normal distribution. 



for r = 1, 2, 3, 4, and Zi is the quantile 
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13. EM Algorithm 

Let the complete-data be Yi, • • • ,Yn with observed values yi, • • • ,yn and the hypothetical 
random variable Zi, - ■ ■ ,Zn. The joint probability density function is such that the marginal 
density of Yi, • • • , is the likelihood of interest. Then, we define a hypothetical complete-data 
distribution for each {Yi,Zi) i = 1, ■ ■ ■ ,n with a joint probability density function in the form 

g{y,z;e) = zajP^y-^-^e-^^yP {1 - e-^^y^Y^-^^, (30) 

where @ = {a, /3,^,9), y > and ^ G N. 

Under the formulation, the E-step of an EM cycle requires the expectation of {Z\Y; Q^^^) where 
Q(^) = (a^^\ (j^^^j^^^O^"^^) is the current estimate (in the rth iteration) of 0. 
The pdf of Z given Y, say g{z\y) is given by 



Thus, its expected value is given by 

5t(l_e-((9!')^)"C"[6»(l-e-('3!')^)° 



E[Z\Y = y] =1 + 



c'\eii-e-(f^y)^)° 



The EM cycle is completed with the M-step by using the maximum likelihood estimation 
over 6, with the missing Z's replaced by their conditional expectations given above. 
The log-likelihood for the complete-data is 

^nivir-- ,yn;zi,--- ,Zn;@) (X J2i=i^og Zi + n[log a + log J + 'y log /3]-Yli=iil3yiV 

+ Er=i(7 - 1) logyi + EILi iog(i - e-(^y^r)^azi - 1) 

-nlogC(e) + Er=i(^i)log^- 
The components of the score function U*{Q) = ( g^, gf-, -gf-) are given by 



da 



= ^+Y:7=i^og{l-e-(M-)zi, 



a/3 ~ /3 iP 2^i=iyi i_e-(/3!'z)T ' 

3^ - ^ +"iog/? -h 2^j=i iogyj - 2^j^i ^_g_(^j,.)7 , 

96> ~ c{e) 61 
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From a nonlinear system of equations U*{Q) = 0, we obtain the iterative procedure of the 
EM algorithm as 

a{t+l) ^ v^n (i) 



^(^+TT + iOg + 2^ .^^ log - ^.^^ ^_^_(^(t)^.)^(*+i) 

where ^(*+^), /3(*+-'^) and are found numerically. Hence, for i = 1, • • • ,n, we have that 



4" = 1 + 



14. Special cases of the EWPS distribution 

In this section we study in detail cases of the EWPS class of distributions. To illustrate 
the flexibility of the distributions, plots of the pdf and hazard function for some values of the 
parameters are presented. 

14- 1- Exponentiated weihull binomial distribution 

The exponentiated weibull binomial distribution is a special case of power series distributions 
with a„ = C^) and C{9) = {9 + 1)™" — 1 {9 > 0), where m (n < m) is the number of replicas. 
Using the cdf in ([7]), the cdf of exponentiated weibull binomial (EWB) distribution is given by 

, (0(l-e-(^2')^)" + l)™-l 

(q(-\ „— (/3j/)T\a I -\\ m—l 

f{y) = ma9jf3^y^-'e-((^yril - ^-(/^.r l^^^ (g + i)i _ 1 ' 

and 

(of-] „— I 1 '\ m—l 

^ ma0T,/37y7-ie-(/3s/r(i _ e-(/3j/)^)"-i ^1^1 f LJLll 

n[y) maujp y e [l e ) ^ -^^m _ (^0^1 _ e-m-')c' + l) m' 

The plots of density and hazard rate function of EWB distribution for some values of a, /?, 7, 
and m = 10 are given in Fig. 2. 

From (I14p , the moment generating function of EWB is 
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Figure 2: Plots of pdf and hazard rate function of EWB for different values a, P,'y , 6 and m = 10. 

A/f ^-^\ — ad ^^^oo v^m v^co {—'^V ( t \inn—l (rn\ (na—l\'"'^(}~^i^ 
MY[t) - 2^j=o Z^n=l Z^j=0 -^{^) ^ [n)\ j /^"^^yT^' 



Q,6)r 1 + - moo / \ / i\ 

14-2. Exponentiated weihull poisson distribution 

The exponentiated weibull poisson distribution is a special case of power series distributions 
with a„ = n!^^ and C{9) = e^ — 1 {9 > Q). Using the cdf in ([7]), the cdf of exponentiated weibull 
poisson (EWP) distribution is given by 

F{y) = ^- 

e^ — 1 

— 1 ' 

and 

Hy) = 



_ ^e{l-e-(M'>)o> 

The plots of density and hazard rate function of EWP distribution for some values of a, /?, 7 
and 9 are given in Fig. 3. 
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Figure 3: Plots of pdf and hazard rate function of EWP for different values a, /3, 7 and 6. 



From p4|) . the moment generating function of EWP is 

M — Y^cx3 v^oo v^oo (— l)i gn-i ^^t^^/na— 1'\ ^(^+~) 

1VIY{1') — ^^e^l) Z^i=0 2^n=l l^j=0 il {n-iyA fS) \ j J 



fe \ 00 00 / ,\ ^ra— 1 



Exponentiated weibull geometric distribution 
The exponentiated weibull geometric distribution is a special case of power series distri- 
butions with a„ = 1 and (7(6*) = e{l - 6)-'^ {0 < 9 < 1). Using the cdf in ([7]), the cdf of 
exponentiated weibull poisson (EWG) distribution is given by 



l_0(l_e-(/32/)^)" ' 
(1 - e)a-f(3^y^-^e-^M'< (1 _ e-i^vry 

~ [I - 9 [I - e-m'')''f ' 

and 

_ (1 - 9)ajf3^y^-'e-(^yy (l - e-(/^s^)")""' 

h{y) - 



[i-e{i- e-(py)'y] [1 - (1 - e-^M-y] • 

The plots of density and hazard rate function of EWG distribution for some values of a, /3, 7 
and 9 are given in Fig. 4. 
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Figure 4: Plots of pdf and hazard rate function of EWG for different values a, /3, 7 and 9. 
From p4|) . the moment generating function of EWP is 

My{t) = a(l - 9) Er=oEn=i Er=o(-iy (T')r U + 1) ^^n9-~Hj + l)^^"'''^- 



-(^+1) 

^y--j y ■ u ^ 

i^.^. Exponentiated weihull logarithmic distribution 

The exponentiated weibull logarithmic distribution is a special case of power series distri- 
butions with a„ = n^^ and C{6) = - log (1 - 6*) (0 < 6* < 1). Using the cdf in ([7]), the cdf of 
exponentiated weibull poisson (EWL) distribution is given by 

log(l-g(l-e-(/^^)"n 
log(l - 6) 

_ Qa^pv^-^e-^f'yr' (1 - e-^f'yr'f-^ 

and 



[^(l-e-(^.)^r-l] log( ,_,(,_^-%.),)J ' 
The plots of density and hazard rate function of EWL distribution for some values of a, /3, 7 and 
Q are given in Fig. 5 



19 



Figure 5: Plots of pdf and hazard rate function of EWL for different values a, P, 7 and 6. 



From p4|) . the moment generating function of EWL is 



Mvit) 



log(i-e) 



ae 



Eoo v^m v^oo / -1 \j 
i=0 2^n=l Z-jj=0\^) 




pi /na~l\ 



ev (1 + f ) 



m 00 / ^ \ 

n=l j=0 \ J / 



& 



in— 1 



log(l-0) 



(i + 1) 




15. Applications of the EWPS distribution 

In this section we present an application of the EWPS to three real data sets. The fit of 
EWG, EWP, and EWL on real data sets is examined by graghical methods using MLEs. They 
are also compared with the EW and Weibull models with respective densities. 
The first data set is given by Barreto-Souza(2009), Morais and Cordeiro on the fatigue life 
(rounded to the nearest thousand cycles) for 67 specimens of Alloy T7987 that failed before 
having accumulated 300 thousand cycles of testing. 

Now, we estimate the parameters of distributions and compare the p-values of Kolmogorov- 
Smirnov test and AIC (Akaike Information Criterion), AD (Anderson-Darling statistic) and 
CM (Cramr-von Mises statistic) for these distributions. 

The empirical scaled TTT transform (Aarset, |11]) and Kaplan-Meier Curve can be used to 
identify the shape of the hazard function. 



20 



Scaled TTT transform 



Kaplan-Meier Curve 




Figure 6: TTT plots and Kaplan-Meier curves of data 1, data 2 and data 3 . 

The TTT plot and Kaplan-Meier curve for the first data in Fig. 6 shows an increasing hazard 
rate function. 

Table 2 lists the MLEs of the parameters, the values of K-S (Kolmogorov-Smirnov) statistic 
with its respective p-value, -21og(L), AIC (Akaike Information Criterion), AD (Anderson-Darling 
statistic) and CM (Cramr-von Mises statistic) for the first data. These values show that the 
EWG, EWL and EW distributions provide a better fit than the EWP and Weibuh for fitting 
the first data. 

We apply the Arderson-Darling (AD) and Cramr-von Mises (CM) statistics, in order to verify 
which distribution fits better to this data. The AD and CM test statistics are described in 
details in Chen and Balakrishnan [3]. In general, the smaller the values of AD and CM, the 
better the fit to the data. As a second application, we consider the data show the fatigue life of 
6061-T6 aluminum coupons cut parallel to the direction of rolling and oscillated at 18 cyclers per 
second. The pooled data, yielding a total of 101 observations, were first analyzed by Birnbaum 
and Saunders (1969). The TTT plot and Kaplan- Meier curve for this data in Fig. 6 shows an 
increasing hazard rate function. 

The MLEs of the parameters, the values of K-S statistic, p-value, -21og(L), AIC, AD and CM 
are listed in Table 3. From these values, we note that the EWG model is better than the EWP, 
EWL, EW and Weibull distributions in terms of fitting to this data. The last data set consists 
101 observations show the stress-rupture life of kevlar 49/epoxy strands which were subjected to 
constant sustained pressure at the 90 stress level until all had failed. The failure times in hours 
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Table 2: MLEs(stds.), K-S statistics, p-values, -21og(L) , AIC, AD and CM for data 1. 



Dist. 


MLEs 




K-S 


value 


-21og(L) 


AIC 


AD 


CM 


EWG 


a = 15.3396 ,/3 
7 = 1.3155 ,6 = 


= 0.0154 
= 0.1860 


0.0486 


0.9974 


695.9917 


703.9917 


0.1968 


0.1029 


EWP 


a. = 20.48 = 
7 = 0.7316 ,e = 


0.0732 
= 13.74 


0.0717 


0.8811 


696.2272 


704.2272 


0.2205 


0.1128 


EWL 


a = 14.0601 ,/3 
7 = 1.3671 ,6 = 


= 0.0158 
= 0.7721 


0.0524 


0.993 


696.8654 


704.8654 


0.2956 


0.1165 


EW 


a = 12.1645 ,/3 
7 = 1.4034 


= 0.0134 


0.0522 


0.9931 


696.0166 


702.0166 


0.19097 


0.1023 


WeibuU 


/3=0.0054 , 7= 


=3.7349 


0.1027 


0.4793 


706.598 


710.598 


1.1684 


0.2541 



Table 3: MLEs(stds.), K-S statistics, p-values, -21og(L) , AIC, AD and CM for data 2. 



Dist. 


MLEs 


K-S 


p- value 


-21og(L) 


AIC 


AD 


CM 


EWG 


a = 8.0516 ,^ = 0.0129 
7 = 2.3695 ,^ = 0.7745 


0.0618 


0.8352 


913.1816 


921.1816 


0.3426 


0.1299 


EWP 


a = 14.022 = 0.0135 
7 = 2.1176 ,^ = 1.059 


0.0791 


0.552 


913.4216 


921.4216 


0.4363 


0.1557 


EWL 


a = 8.9561 ,^ = 0.01143 
7 = 2.4247 , e = 0.2769 


0.0832 


0.4867 


913.7988 


921.7988 


0.5413 


0.1729 


EW 


a = 8.072 ,/3 = 0.0108 
7 = 2.5872 


0.082 


0.5049 


913.498 


919.498 


0.4597 


0.1616 


WeibuU 


/3=0.0069 , 7=6.0347 


0.1234 


0.0923 


926.9108 


930.9108 


1.755 


0.3657 
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Figure 7: Fitted cdf and survival function of the EWG, EWP, EWL, EW and Weibull distributions for the data 
sets corresponding to Table 2. 

are shown in Andrews and Herzberg and Barlow et al. |4]. The TTT plot and Kaplan-Meier 
curve for this data in Fig. 6 shows bathtub-shaped hazard rate function. The MLEs of the 
parameters, the values of K-S statistic, p-value, -21og(L), AIC, AD and CM are listed in Table 
4. From these values, we note that the EWG and EWP models are better than the EW and 
Weibull distributions in terms of fitting to this data. 



Table 4: MLEs(stds.), K-S statistics, p-values, -21og(L) , AIC, AD and CM for data 3. 



Dist. 


MLEs 




K-S 


p- value 


-21og(L) 


AIC 


AD 


CM 


EWG 


a = 1.0921 ,/3 

7 = 0.661 ,e = 


^ 3.1202 
0.7559 


0.0724 


0.6657 


203.66 


211.66 


0.7842 


0.2019 


EWP 


a = 0.8589 ,/3 
7 = 0.8717 ,6 = 


= 1.3032 
= 1.2661 


0.0725 


0.6638 


204.6174 


212.6174 


0.8409 


0.2182 


EWL 


a ^ 2.4513 ,/3 
7 = 0.4978 , e = 


= 17.0129 
= 0.9918 


0.0898 


0.3893 


202.4622 


210.4622 


0.8643 


0.2455 


EW 


a = 0.7929 ,/3 
7 = 1.0604 


= 0.8210 


0.0844 


0.468 


205.5743 


211.5743 


0.9554 


0.2473 


Weibull 


;3=1.0101 , 7 


=0.9259 


0.0906 


0.3778 


205.9536 


209.9536 


1.1221 


0.2789 



Plots of the densities and cumulative distribution functions of the EWG, EWP, EWL, EW 
and Weibull models fitted to the data sets corresponding to Tables 2, 3 and 4, respectively, are 
given in Fig. 7, 8 and 9. 
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Figure 8: Fitted cdf and survival function of the EWG, EWP, EWL, EW and Weibull distributions for the data 
sets corresponding to Table 3. 




Figure 9: Fitted cdf and survival function of the EWG, EWP, EWL, EW and Weibull distributions for the data 
sets corresponding to Table 4. 
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16. Conclusion 



Appendix 

The elements of the 4x4 observed information matrix 7„ (G) are given by 



J _ _ n I Y^n 6>(log(l-e-(/3i^i)^))^(l-e-(^^'i)^)"C"(6>(l-e-(^t'i)^)") 
^aa — ^ -t- l^i=l C"(6>(l-e-(''!'<)'^)«) 

6l2[(l-e-(^!'i)^)" log(l-e-('3!'i)^)]2C""(6l(l-e-W!'i)^)«) 
+ Z^i=l C'(6»(l-e-('3!'»)^)'») 

Z^i=l (C"(6»(l-e-(^2'i)'')«))2 ' 



-'a/? — Z^i=l l_e-(/3yj)7 

6>7/3T-ii/iTe-tfi^»)^(l-e-(^i'i)^)"-iC"(9(l-e-(^^»)^)") 
+ 2^1=1 c"(6»(l-e-('32'i)^)") 

+ 2^1=1 C"(6>(l-e-(^2'i)'')") 

2^i=l C"(6l(l-e-('3!'i)''')«) 

Z^i=l (C"(6»(l-e-(^2'i)'^)°))2 



2 





/3Tj/iTlog(/3j/i)e-('3s'i)'' 




l_e-(/3yi)'V 




+Er=i 


e/9Tj/i'i'log(/3j/i)e-(''2'i)^(l-e-('3!'i)^)"-iC"(e(l-e-(''3'<)'')«) 




C'(6l(l-e-('9!'i)^)") 






0a/3TyiTlog(/3i/i)e-(''w)^log(l-e-W2'»)^)(l-e-(''w)'')"-iC"(6»(l-e- 


-(/3!/i)^)«) 


C"(6>(l-e-('5!'<)^)«) 






e2a/3T2/iTlog(/9j/i)e-(''!'i)^log(l-e-(''!'i)^)(l-e-(^!'i)^)2«-iC""(6l(l- 


-e-Wfi)'')™) 


C'{6»(l-e-(''!'i)^)") 




2^1=1 


6»2a/3Tyi^log(/3yi)e-(''2'i)^ log(l-e-(''2'i)^)(l-e-('3!'«r' )2<^-i (C"(6I(1- 


_e-(/5ai)^)"))' 


(C"(6»(l-e-(/'!'«)^)«))2 





^ _ log(l-e-(^a^')^)(l-e-tf^'i)^)"C'"(e(l-e-(^i^»)'^>) 
— 2^1=1 C''(g(i_e-(^yi)'^)a) 

v^n eiog(l-e-(^i'')"^)(l-e-(^?'i)^)2"C'"(6i(l-e-(^^i)^)°) 

+ 2^j=i c'(e(i-e-('3!'i)^)") 

eiog(l-e-(^!^i)^)(l-e-tf!^i)]')2"_(C"(6'(l-e-(^^<»)^)"))2 
2^1=1 (C"(6»(l-e-(^!'i''^)«))2 
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-'/S/S ~ ^ Z^i=l l\I ^JP yi ^ Z^i=l i-e-iPyir 

2^i=l i_e-(/3ai)^ 2^j=l (i_.e-('9!'i)^)2 

, v^n ea7(7-l)i/iT/3T-^e-(^^i)^(l-e-(^i^»)2)°-^C"(e(l-e-tf^i)^)") 
2^1=1 c'/(6)(i_e-(/3!/i)''')«) 

_ 6'a72i/,-\.J- -(. (7" (^(l-(. 

2^1=1 c'{e{i-e~'^i^yi>'')") 
+ 2^1=1 

En 
1=1 



C'(6>(l-e 


-(/3!/i)-')<:<) 




e2Q;272yi27/327-2 (g-(/3j,i)T )2 (i_g- 


-(^!,j)T)2a-2(^///(g(l_g-0 




c'(e»(i-e- 


(/3ai)T)«) 




e2a2722,,27^27-2(g-(,3y,)T)2(i_g- 


-('3!'i)'^)2"-2(C"(6l(l-e-( 


/5j'i)'')«))2 



I/Sj = ^ - ELi 7/3^- log (PVi) - Er=i 

"I" 2^1=1 i-e-C/Sai)"' "T 2^i=l i_e-(/3»i)"' 

_ v-n (a-l)7^^^-^yi^^log(/3yi)e-(^^»)"' _ ,r-n (0-1)7^27-1^.27 ipg (^^ ^j)T)2 
2^1=1 i-e-C/SfiP 2^i=l (i_e-(/3!/j)T)2 

, -^n ea7i/iT/3T-Mog(/3yi)e-(^^»)^(l-e-(^i'»)^)"-^C"(6'(l-e-(^^»)^)°) 

, -^n 6>ai/jT/3T-ie-tf^»)^(l-e-(<^^»)^)"-iC"(e(l-e-tf^»)^)") 

En daryi-'J-' Mog(.ji/,)c ^'y.'^l-c ^C;"((9(l-c '"'.''')") 

1=1 c'/(0(i_e-(/3ai)''')<^) 

, >pn (a-l)6'a7yi 27/327-1 j;i)T)2(^ -(/3^j)T)a-2(g//(g(i_g-tf3;i)T)a-) 

2^1=1 c'/(6)(i_e-(/3!'iF)a) 

, v^n e^a^7yi^T/32T--ilog(/3^,)(e-(/i'!'.)"')2(l-e-tf^^»)^)2°-2c"'(e(l-e-tf^^i)^)°) 
2^1=1 C'/(e(i_e-(/3yi)T)a) 

_-^n 6>2a27j/i27/327-lipg(^^.)(g-(/3yi)T)2(^_e-(/3^,)^-)2a-2(g,»(g(i_g-(j33,j)T)«^)2 

2^1=1 (c/(6)(i_e-(/3ai)'')"))2 ' 



J _ v^n a7yiT/3T-^e-(^!^i)^(l-e-tf^i)"^)"-iC"(6>(l-e-(^^»)'^)") 
'/3e - 2^i=i C"(e(l-e-('^!'i) ')") 

, v^n ea7^iT^T-ie-(^^»)^(l-e-(<^^»)^)2"-iC'"(6>(l-e-(<^^»)^)") 
2^1=1 (^/(^(i-e-C/Sj/iP)") 

_ v^n 6'a7yiT/3T-^e-tf^')^(l-e-(^i^»)^)2"-i(C"(6)(l-e-(^i'»)^)"))^ 
2^1=1 (C'(6'(l-e-(''2'i)'')'»))2 
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(a-l)/3^Tyi^^( log (/3y,))^e-(g^»)^ y^n (a-l)/3^^i/i^T( log (/3i/i))^e-(^i'i)'"' 
Z^i=l i_e-(/3ai)^ 2^1=1 (i_e- )2 

6>a/3TyiT(log(/3yi))^e-tf^i)^(l-e-tf^i)^)°-iC"(e(l-e-tf^i)^)°) 



C"(6l(l-e-('3!'j)^')«) 



1=1 r;'rfli'i-fi-('3!/i) 



, v;^n 6ia(a-l)/3^Ti/i^T(log (/3yi))^(e-(^^i)^)^(l-e-tf^»)^)°-^C"(6>(l-e-(^^i)^)«) 
+ 2^1=1 



C'{e{l-e- 










-2C""(6>(l-e-(''!'i)^)«) 


C'(e(i-e-' 


:/5j/i)T)c«) 






.e-(/3!'i)^)2«- 


-2(C"(6»(l-e-W!/i)'')«))2 



En ^ 
i=l (C'(6'(l-e-('32'i''')«))2 



r _ v^n a/3TyjTlog(/3yi)e-(^i^»)^(l-e-(^s'»)'^)"-iC"(6>(l-e-(g^i)^)") 
7f - Z^i=l C"(e(l-e-(^!'i)'')°) 

, v^n 6>a/3T'j/iT'log(/3j/i)e-tf!'i)^(l-e-(^^^i)^)^"-^C"'(6>(l-e-Cay,)' )a) 
+ 2^1=1 C"(6»(l-e-(''!'<)''')«) 

_-^n 6>a/3TyiTlog (/3yi)e-(^i'»)^(l-e-(^i^»)^)2°-i(C"(6>(l-e-(^«i)^)°))^ 



_ _ n I v^n (l~e-<-l^yi'l^)^''C"'{e{l-e-(f^yO'^)°') 

J-ee — -ffi -r i^i=i (;/(e(i_e-(/3Hi)'')c) 

_ (l-e-(^^'i)^)^"(C"(e(l-e-<^^')^)°))^ _ nC"(6>) , n(C'(6l))2 

z^i=i (c"(e(i-e-('3!'i'^)«))2 c{e) (C(e))^ 
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